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We present a model for wave propagation in a monolayer of spheres on an elastic substrate. The 
model, which considers sagittally polarized waves, includes: horizontal, vertical, and rotational de¬ 
grees of freedom; normal and shear coupling between the spheres and substrate, as well as between 
adjacent spheres; and the effects of wave propagation in the elastic substrate. For a monolayer of 
interacting spheres, we find three contact resonances, whose frequencies are given by simple closed- 
form expressions. For a monolayer of isolated spheres, only two resonances are present. The contact 
resonances couple to surface acoustic waves in the substrate, leading to mode hybridization and 
“avoided crossing” phenomena. We present dispersion curves for a monolayer of silica microspheres 
on a silica substrate, assuming adhesive, Hertzian interactions, and compare calculations using an 
effective medium approximation to a discrete model of a monolayer on a rigid substrate. While 
the effective medium model does not account for discrete lattice effects at short wavelengths, we 
find that it is well suited for describing the interaction between the monolayer and substrate in the 
long wavelength limit. We suggest that a complete picture of the dynamics of a discrete mono- 
layer adhered to an elastic substrate can be found using a combination of the results presented 
for the discrete and effective medium descriptions. This model is potentially scalable for use with 
both micro- and macroscale systems, and offers the prospect of experimentally extracting contact 
stiffnesses from measurements of acoustic dispersion. 

PACS numbers: 45.70.-n, 62.30.-l-d, 63.20.D-, 68.35.Np 


I. INTRODUCTION 

Granular media are simultaneously one of the most 
common and complex forms of matter on Earth. This 
complexity stems, in part, from heterogeneous structure 
and highly nonlinear particulate interactions m- Over 
the past thirty years, mechanical wave propagation in 
ordered granular media has become an active field of 
research as it provides a setting for the broader under¬ 
standing of granular media dynamics [3]. Ordered gran¬ 
ular media have also been shown to enable a wide ar¬ 
ray of novel passive wave tailoring devices that leverage 
the nonlinear response stemming from the Hertzian re¬ 
lationship between elastic particles in contact [3 IS], in 
conjunction with dispersion induced by periodicity [5] or 
local resonances [7|. 

Experimental configurations used to study mechanical 
wave propagation in ordered granular media typically in¬ 
volve spherical particles confined by elastic media. This 
type of arrangement is particularly common in one- and 
two-dimensional configurations, and includes macro- to 
microscale particles. For example, at the macroscale, 
elastic rod structures, tracks, and tubes have been used to 
confine the particles in one-dimensional and quasi- 
one-dimensional uni configurations, and elastic plates 
have been used in two dimensions HD. More recently, 
the dynamics of a two-dimensional monolayer of 1 ^m 
silica particles adhered to an elastic substrate was stud¬ 
ied using a laser ultrasonic technique [12) . 

Analytical models used to describe the dynamics of 
these systems typically only include the interaction be¬ 
tween the particles (often just the normal Hertzian con¬ 
tact interaction) and disregard the effect of the substrate. 


In reality, even for the simple case of a particle monolayer 
on a substrate, more complex dynamics involving inter¬ 
actions between the particles and elastic waves in the 
substrate should be expected. Indeed, a recent experi¬ 
ment m showed that a monolayer of microspheres ad¬ 
hered to a substrate strongly interacts with Rayleigh sur¬ 
face waves in the substrate, leading to the hybridization 
between Rayleigh waves and a microsphere contact reso¬ 
nance. The results of this experiment were analyzed with 
a simple model involving only vertical (normal to the 
substrate surface) vibrations of isolated particles, follow¬ 
ing the approach adopted in earlier theoretical works on 
the interaction of surface oscillators with Rayleigh waves 
[nun]. However, in reality, the particle motion is not 
confined to the vertical direction, and the Rayleigh wave 
has a significant horizontal component. Furthermore, the 
interaction between neighboring particles in contact is 
expected to significantly influence the dynamics. 

A notable theoretical work [15] provided a model for 
the dynamics of a monolayer adhered to an elastic sub¬ 
strate which accounted for both normal and horizontal 
motion and interaction between the particles. However, 
this study did not take into account particle rotation. A 
more recent study m demonstrated that the rotational 
degree of freedom has a profound effect on the dynamics 
of granular monolayers. However, the analysis of mono- 
layers on substrates in Ref. m only accounted for nor¬ 
mal contact interactions between the particles and the 
substrate, and the substrate was considered rigid. 

The aim of this work is to provide a theoretical model 
for the contact-based dynamics of a two-dimensional 
layer of spheres on a substrate, accounting for the elas¬ 
ticity of the substrate, translational and rotational mo- 


2 


tion of the spheres, and both normal and shear stiffnesses 
of sphere-to-sphere and sphere-substrate contacts. We 
focus on a system with microscale particles that inter¬ 
act with each other and with the substrate via van der 
Waals adhesion forces. Rather than postulate the con¬ 
tact stiffness constants, we derive them from Hertzian 
contact models. This imposes certain constraints on the 
values of the constants: for example, the ratio of the nor¬ 
mal and shear contact stiffness between the spheres is a 
constant only weakly dependent on Poissons ratio. We 
consider contact-based modes having frequencies signif¬ 
icantly below the intrinsic spheroidal vibrational modes 
of the spheres, such that they can be described as spring- 
mass oscillators. Furthermore, we focus on dynamics in¬ 
volving particle and substrate displacements in the sagit¬ 
tal plane, as would be detectable in a laser-based exper¬ 
iment, such as that of Ref. [12] . 

We start with the case of a rigid substrate, where we 
find three eigenmodes involving vertical, horizontal and 
rotational motion of the spheres. In the long-wavelength 
limit these modes yield three contact resonances, for 
which simple analytical expressions are obtained. One 
of the resonances only involves motion of the spheres 
normal to the substrate surface, whereas the other two 
involve mixed horizontal-rotational motion. We then 
present our effective medium model, which describes the 
interaction between the spheres and the substrate. The 
results show that the contact resonances interact with 
Rayleigh surface waves, which leads to mode hybridiza¬ 
tion and avoided crossings. We discuss the behavior for 
cases involving both isolated (non-touching) and inter¬ 
acting spheres, and demonstrate the important role of 
rotations in both cases. We also examine the validity of 
the effective medium approximation, by comparing the 
calculations using discrete and effective medium models. 
Finally, we discuss the implications of our findings for 
past and future studies on granular monolayer systems. 


II. MODEL 

We consider a monolayer of elastic spheres on a sub¬ 
strate, which can be either close-packed and in contact, or 
isolated, as shown in Fig. [^a). In either case, the spheres 
are assumed to form a square lattice, with the wave prop¬ 
agation direction aligned with the lattice vector, as shown 
in Fig. 0 b)- We model the layer as an infinite lattice 
of rigid spheres with diameter D = 2R and mass to, 
coupled to a semi-infinite, isotropic elastic substrate by 
normal and shear stiffnesses RTat and RTg, and to nearest- 
neighbor spheres by stiffnesses Gm and G 5 , as schemat¬ 
ically shown in Fig. 0c). The subscript N corresponds 
to forces acting normal to the surface of the sphere, and 
S to forces acting transverse to the surface of the sphere. 
The shear springs generate an associated torque about 
the sphere center, while the normal springs do not. The 
absolute horizontal, vertical, and angular displacements 



FIG. 1: (a) Side-view schematic of an amplified wave profile 
for isolated and interacting spheres, (b) Top-down view of 
the square-packed monolayer, with the arrow indicating the 
direction of wave propagation, (c) Schematic for the model 
of a monolayer of spheres coupled to an elastic halfspace. 

of sphere j from its equilibrium state are given by Qj, 
Zj, and 9j, respectively, and the displacements of the 
substrate are denoted by u(x,z), corresponding to dis¬ 
placement in the x-direction, and w{x, z), corresponding 
to displacements in the z-direction. 

A. Contact Stiffness 

We derive the stiffnesses iFjv, Ks, Gjq, and Gg using 
Derjaguin-Muller-Toporov (DMT) [111113 and Mindlin 
contact models [21]. The DMT theory is typically ap¬ 
plicable in weakly-adhesive systems with small, stiff par¬ 
ticles |Hj, and assumes that the deformation profile is 
Hertzian. The Mindlin model describes the shear stiffness 
of particles in contact, assuming an applied normal force 
[21j . At the microscale, adhesive contact models have 
been explored experimentally in the quasi-static regime 
using atomic force microscopy and nanoindentation ap¬ 
proaches [ 53 . 

For contact between two spheres (or a sphere and a 
halfspace) having elastic moduli Ex and ^ 2 , and Pois¬ 
son’s ratios vi and 1 ^ 2 , the Hertzian restoring elastic force 
Fjq corresponding to displacement 6n of the particle cen¬ 
ter in the direction normal to the contact surface is given 

by 

Ft, = ( 1 ) 

where Rc is the effective radius of contact (equal to R 
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for sphere-sphere contacts and R/2 for sphere-halfspace 
contacts), and E* = [(1 — v1)/Ei -I- (1 — is the 

effective modulus. Considering the DMT adhesive force 
Edmt = 2ttwRc acting normal to the contact surface 
[ini [ 20 ] (where w is the work of adhesion between two 
surfaces), the net normal force is given by 

FN,net = Fn — Edmt- ( 2 ) 


To describe the shear contact, we utilize the Mindlin 
model ED, which assumes small relative displacements 
and no slip at the contact surface. For two elastic bod¬ 
ies with shear moduli Gi and G 2 , the restoring elastic 
force Fs to displacement Ss of the particle center in the 
direction transverse to the contact normal is given by | 21 ] 

Fs = 8G*Rl^^Ss6^J\ (3) 


where G* = [(2 — j^i)/Gi-I-(2 — z/ 2 )/G 2 ] ^ is the effective 
shear modulus. Here, the factor of arises from the 
Hertzian relation between the contact radius and Fj^. 

By substituting the relative displacements 6n = Z—wq 
and Ss = Q — uq + R9 into Eqs. and and lin¬ 
earizing about the equilibrium configuration of S^ q = 

[ 3 Fi 5 Mr/( 4 F;*i?y")] 2/3 and Js.o = 0 , we derive linearized 
normal and shear contact stiffnesses 


Km = 6E* RrF, 


,1/3 


cJ^^DMT , 


/ Q Q* \ 

^S = 8y-—RcFDMTj 


1/3 


(4) 


with Gjm and G 5 given by equations of the same form, 
but with F*, G*, Rc, and Fjjmt adjusted for sphere- 
sphere contacts. In the special case where the spheres 
and substrate are composed of the same material, the 
relative magnitudes of the stiffness constants are deter¬ 
mined exclusively by Poisson’s ratio v of the material, 

Gn = 

Ks = u* Km (5) 

Gs = 2 -'^I^v*Km, 

where v* = 2(1 - jz)/(2 - v). 


B. Equations of Motion of the Spheres 


Assuming small displacements (i.e. Qj, Zj^ and RQ^ 
are much less than D), the sphere obeys the equations 


of motion 

ttiQ j — Fs(^Qj 4” ^^j') 

+ GM{Qj+i — 2Qj + Qj-i) 

mZj = — KM{Zj — wqj) 

+ Gs[Zj+i — 2Zj + Zj_i — i?(0j_|_i — dj_i)] 

16 j = — KsR{Qj — uoj + R6j) 

— GsR\R{6jj^i -\- 26j 6j—\) 

- {Zj+i - Zj_i)\, 

( 6 ) 


where mqj and wqj are horizontal and ver¬ 
tical displacements of the substrate sur¬ 
face at the point of contact, respectively. 


C. Effective Medium Approximation 


Considering wavelengths much longer than the sphere 
diameter, we treat the monolayer as an effective contin¬ 
uous medium. By substituting the center difference for¬ 
mulas [(Oi-ei - (•)j-i]/(2-D) « d{-)/dx and - 

2(- )j + (• )j-i]/{D^) « d'^i- )/dx^ into Eq. ([^, we arrive 
at the equations of motion of the monolayer in effective 
medium form: 


m 


d^Q 

dt^ 


d^Z 


d^e 


-Ks{Q-uo + R6) 

+ AGnR^^ 

OX^ 

- Kn{Z - Wo) 

,d^Z 36 

+ 4G5i? (^-^) 

- KsR{Q -U 0 + R6) 

- 4Gsi?^(i ?^—3 + 6+ —). 


dZ^ 


dx^ 


dx ' 


(7) 


The coupling between the monolayer and substrate is de¬ 
scribed by the following boundary conditions at the sur¬ 
face z = 0 , which describe the average force acting on 
the surface due to the motion of the spheres: 


'X zx 

O’zz 


A 


(Q — Ug -h R6) 


A 


(Z - Wo), 


( 8 ) 


where azx and a^z are components of the elastic 
stress tensor m and A = is the area of a 

primitive unit cell in our square-packed monolayer. 
The combination of Eq. Q and the linear elastic 
wave equations describing waves in the substrate 
lljj . coupled by the boundary conditions of Eq. 
®, comprises the complete effective medium model. 
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III. DISPERSION RELATIONS 
A. Rigid Substrate 

1. Discrete Model 

We first consider the discrete model, which accurately 
captures the structural periodicity of the monolayer. We 
assume spatially-discrete traveling wave solutions of the 
form (with similar terms for the other dis¬ 

placements) and set the displacements of the substrate 


surface uqj and wqj to zero (rigid substrate). Here, (•) 
is the amplitude of a plane wave in the displacement vari¬ 
able (•), w is the angular frequency, and k is the wave 
number. After algebraic manipulation, Eq. ([^ is re¬ 
duced to a homogeneous system of three linear algebraic 
equations in the amplitudes (•). This system possesses 
non-trivial solutions only for pairs of k and w that cause 
the determinant of the system to vanish. Enforcing this 
condition, we arrive at the following dispersion relation, 
where the three rows of the determinant correspond to 
the three equations of Eq. (§: 


^(1 - cos{kD)) + (j)sujl 
0 

Ru^l 


0 


RujI 


^(1 - cos(fc£)))-I--|^isin(fcD) 
^is\n{kD) ^[^(l-cos(fcD))H 




= 0 , 


(9) 


where (f)N = 1 — (f)s = 1 — (j)g = I — 

Cm = yjGN/tn{2R) and cs = ^JGs/m{2R) are the 
longitudinal and transverse long-wavelength sound 
speeds of the discrete monolayer, Cg = —ruR’^Cg/I, 
ujjf = Km/tti, ujg = Ks/m, and ujg = {Ks + 4:Gs)R^/I- 


spatially-continuous traveling wave solutions of the form 
Q^i{ujt-kx) similar terms for the other displace¬ 

ments) into Q, Z, and 0 in Eq. Q with uq = wq = 0. 
Following the same procedure as in the discrete model, we 
arrive at the dispersion relation for the effective medium 
model: 


2. Effective Medium 

In the effective medium model, which approximates 
the discrete system at long wavelengths, we substitute 


c%k'^ + (j)S^s 0 

0 c%k'^ + 4>NU}l! 

Rujg ikcg 


RujI 

—ikcg 

^{c^k^ + (l)gU}l) 


= 0 . 


( 10 ) 


It is particularly instructive to examine the behavior of 
the effective medium model in the long wavelength limit 
A: —>■ 0. In this limit, the displacements vary slowly in 
space, and the spatial derivative terms of Eq. Q may 
be neglected. For the case of a rigid base, Eq. Q then 
reduces to the form 


The equation for Z decouples from the other two equa¬ 
tions and yields a vertical vibrational mode. The two 
other equations remain coupled, yielding two modes con¬ 
taining both horizontal and rotational motion. Using the 
moment of inertia of a solid sphere I = (2/5)mi?^, we 


d^Q 
™ dt^ 
d^Z 


d^e 


Ks{Q + Re) 


KmZ 

KsR{Q + R6) - AGsR^e. 


( 11 ) 
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find three resonance frequencies 
1/2 



Kn 

UJn = 

m 


\fKs 

^RH = 

\4m 


IfKg 

^HR = 

\4m 


49 


(2O7 + 7 + V4OO72 + I2O7 
(2O7 + 7 - V4OO72 + I2O7 + 49) 


1/2 

1/2 

( 12 ) 


where 7 = Gs/Kg. Here, lum corresponds to a mode 
with exclusively vertical motion, described by Eq. 0 - 
The other two modes lurh and ujhr exhibit both ro¬ 
tational and horizontal (but not vertical) motion, with 
relative amplitudes determined by 7 . The higher of the 
two horizontal-rotational modes is predominantly rota¬ 
tional and the lower is predominantly horizontal, hence 
we have used the notations lorh and ujhr, where the 
first letter in the subscript denotes the dominant mo¬ 
tion. If the spheres and substrate are made of the 
same material, then, by using Eq. ([^, we can re¬ 
late the horizontal-rotational frequencies of Eq. (12) 


to the vertical resonance frequency, with the expressions 
^RH = 3.018;/* ^ Wat and ojrr = 0.832;/* ' ujn. 

In the limiting case of isolated spheres (described by 
7 = 0), ujRH and uirr of Eq. (fl^ become ojrhjso = 


yJY/2ws and ojhrjso = 0, respectively. For identical 
materials, ujrh.Iso = \/7;/*/2wjv- The dependence of 
ujRH and ujhh on 7 is shown in Fig. (a), where it can 
be seen that ujrr originates at ojrhjso for 7 = 0 and 
grows unbounded, while lurr originates at ujhr, iso = 0 , 
and approaches uig asymptotically. In Fig. [^(b), we plot 
the horizontal and rotational displacement amplitudes as 
functions of 7 for these two modes. Different signs of 
the rotational amplitude indicate that the lurh and lurr 
modes have different motion patterns. In the former, a 
positive displacement is accompanied by a counterclock¬ 
wise rotation, while in the latter, it is accompanied by a 
clockwise rotation. 

We note that the zero-frequency mode, ujrrjso^ corre¬ 
sponds to the rolling motion of an isolated sphere. With 
the inclusion of a bending rigidity, the sphere would not 
be allowed to freely roll, and instead would undergo rock¬ 
ing motion of a finite frequency. While non-zero bending 
rigidity is expected to exist in real systems (for instance, 
a microsphere adhered to a substrate does not freely roll), 
the frequency of resulting rocking vibrations is predicted 
to be orders of magnitude lower than the other contact 
resonances discussed here |24j . Bending rigidity would 
thus act as a small perturbation to the predictions of our 
model, and we do not include it in our analysis. 

To illustrate the importance of particle rotations in 
the model, we note that in the limiting case of / —>■ 00 , 
when there is no rotation, Eq. 0 yields two resonances: 
a vertical resonance with frequency lor, and a horizon¬ 
tal resonance having frequency tog. For isolated spheres, 
the effect of rotations can be thought of as a reduction 




FIG. 2: (a) Resonance frequencies uirh (red line) and ujhr 
(black line) as functions of the stiffness ratio 7 . (b) Displace¬ 
ment amplitudes of the resonant modes with frequencies ojrh 
(red lines) and ujhr (black lines), as functions of the stiffness 
ratio 7 . Solid and dotted lines correspond to Q and R6, re¬ 
spectively. For each resonance, the amplitudes are normalized 
such that the sum of squares is unity. The positive sign of R9 
corresponds to counterclockwise rotation. 


of the “effective mass” of the sphere to {2(7)m, which 
increases the horizontal resonance frequency. For inter¬ 
acting spheres, on the other hand, rotations drastically 
change the dynamics, yielding two horizontal-rotational 
modes whose freqnencies depend on the relative strengths 
of the sphere-to-sphere and sphere-substrate interactions. 


B. Elastic Substrate 

As in the case of the effective medium approximation 
for a rigid substrate, we assume traveling wave solutions 
of the form Q^^Gt-kx) similar terms for the other 

displacements) into Q, Z, and 9 in Eq. Q. Likewise, we 
express the variables Uq, Wq, a^x <^zz in terms of surface 
wave solutions for the elastic potentials [T7] cj){x,z,t) = 

^^kaz+i{uR-kx) = ^e'^P^+^Gt-kx) ^ 

substitute these expressions into Eq. Q and Eq. 

Here, (•) is the amplitude of a plane wave in the displace¬ 
ment or potential variable (•), a = ^ 

P = VT — uj'^/, and cr and ct are the longitu¬ 
dinal and transverse sound speeds of the substrate, re¬ 
spectively. After algebraic manipulation, Eq. 0 and 
Eq. ( ^ are reduced to a homogeneous system of five lin¬ 
ear algebraic equations in the five plane wave amplitudes 
(•), with coefficients depending on k and to. We reach 
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the dispersion relation by seeking nontrivial solutions of cause the determinant of the following coefficient matrix 

this system, which exist only for pairs of k and oj that to vanish: 


ikujg 

kf3ujg 


0 

—kaujff 

ikojjf 

0 

Cgk'^ + (j)N0J% 

ikRujg 

kRPujg 

R^s 

ikcg 

1 + ^2 

-2i/3 

0 


—2ia 


+ wl) 

0 


m V r 


R^s 

—ikcg 


0^ 

— m 


pAc^k^ *^*^1 
0 


= 0 , 


(13) 


where p is the density of the substrate, and A = is the 
area of a primitive unit cell in our square-packed mono- 
layer. We note that the coupling between the spheres and 
the substrate is represented by elements (4,4), (4, 5), and 
(5, 3) of the matrix in Eq. Thus, the strength of the 
coupling can be quantified by the ratio m/{pA); if this 
term is made to vanish (e.g. by making the mass of each 
sphere much less than that of the portion of the substrate 
below it, extending to the depth of material influenced by 


Rayleigh waves), then the substrate and monolayer will 
be effectively decoupled. We note that if rotations are 
disregarded (e.g. by letting I oo), Eq. (13) reduces to 
the same form as that of the adsorbed monolayer of Ref. 

m- 


It is instructive to consider the long-wave limit when 
the spatial derivatives in Eq. ([^ can be disregarded. In 
this case, we find the simplified dispersion relation 


ikuig 

k/3ujg 

•PsUJs 

0 


—kaujj^ 

ikuj'lf 

0 


0 

ikRijjg 

kRjSujg 

Rul 

0 


1 + 02 

-2i0 

0 

^n) 

0 

—2ia 

-(1 + 0") 


0 

0 


Eor isolated spheres, there is no approximation in Eq. 


(14) with respect to Eq. (13), because in this case the 
terms generated by the spatial derivatives in Eq. ([^ are 
identically zero. Eor interacting spheres, the accuracy 
of dispersion relations calculated with Eq. ( [l4| ) will be 
assessed below by a comparison with results obtained 
with Eq. (13). We will see that Eq. (14) essentially 
describes the interaction of contact resonances given by 
Eq. (11) with Rayleigh surface waves. 


IV. NUMERICAL RESULTS AND DISCUSSION 


In the following calculations, we consider silica spheres 
of 1.08 /im diameter on a silica substrate, and use the 
elastic constants (Ref. [l^) E = 73 GPa, = 0.17, and 
work of adhesion (Ref. [20]) w = 0.063 J/m^. 


A. Rigid Substrate 

We plot numerical solutions of Eq. § , to obtain the 
dispersion curves for the discrete model of interacting 
spheres on a rigid base, as shown in Fig. ^a). In our 
description of a rigid substrate, we assume that no elastic 
waves propagate in the substrate, but allow local defor¬ 
mation at the points of contact for the purpose of the 
contact stiffness calculation; this preserves the same con¬ 
tact stiffnesses as in the elastic substrate analysis. We 
note that due to the periodicity of the system, all three 
branches have zero-group velocities at the edge of the 
first irreducible Brillouin zone |18j of the monolayer. 

By substituting the solutions shown in Fig. da) into 
the coefficient matrix of the corresponding algebraic sys¬ 
tem, we numerically determine the amplitudes of the 
sphere displacements, which we plot in Fig. ib-d). By 
comparing the calculated displacements of with the dis¬ 
persion curves, we see that each branch takes on the 
character of its respective contact resonance in the limit 
k —>■ 0. One can see that each of the three modes gener¬ 
ally involves both vertical and horizontal, as well as rota- 
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FIG. 3: (a) Dispersion relation of a discrete monolayer ad¬ 
hered to a rigid base. Blue solid and red dotted lines de¬ 
note, respectively, discrete and effective medium monolayers. 
Black dashed lines denote the contact resonances, (b)-(d) rel¬ 
ative amplitudes of the displacement variables Q (black dot¬ 
ted lines), Z (red dotted lines), and R9 (blue dotted lines), 
corresponding to the branches of the same numeral for the 
dispersion of the discrete monolayer adhered to the rigid base 
shown in (a). The amplitudes are normalized such that the 
sum of the squares is unity. 


tional motion (albeit the rotational component of mode 
II is quite small). The existence of the three modes with 
mixed displacements is a consequence of the inclusion of 
the rotational degree of freedom: without rotations, there 
would be two modes, one vertical and one horizontal. 

We note that, in the special case Ks = 0, the mode 
originating at ujhr becomes purely horizontal and decou¬ 
ples from the other two modes. The remaining modes 
(characterized by vertical translation and rotation) are 
generally consistent with the results of Ref. [T^, for 
the case of normal contact with a rigid surface and no 
bending rigidity. Since Ref. m considered hexagonal 
packing, the behavior is analogous at long wavelengths, 
but diverges at short wavelengths due to discrete lattice 
effects. 

The dotted lines in Fig. ia) show dispersion curves 



FIG. 4: Dispersion relation of a monolayer adhered to a rigid 
base, using an extended plotting range. Blue solid, red dotted, 
and black dashed lines are the same as in Fig. I^a). Black 
dash-dotted lines denote the long-wavelength longitudinal and 
transverse sounds speeds of the monolayer. 


calculated with the effective medium model as per Eq. 
(10). The effective medium approximation yields accu¬ 
rate results at long wavelengths but fails at shorter wave¬ 
lengths with the unphysical behavior of the first mode, 
whose frequency goes to zero. At even shorter wave¬ 
lengths, as shown in Fig. the effective medium disper¬ 
sion curves of modes II and III asymptotically approach 
straight lines with slopes given by the longitudinal and 
transverse sound speeds in the monolayer. This asymp¬ 
totic behavior has been described by Kosevich and Syrkin 
m- However, as can be seen from the dispersion curves 
generated using the discrete model in Fig. this asymp¬ 
totic behavior does not occur in our system due to the 
spatial periodicity of the monolayer. As a result, the in¬ 
clusion of the first- and second-order spatial derivative 
terms of Eq. ([^ does not deliver much additional under¬ 
standing of the dynamics of our system. 


B. Elastic Substrate 


1. Isolated Spheres 


We numerically solve Eq. (14) for the isolated spheres 
case using Gs = 0 and all other parameters derived 
in Sec. IIA and plot the resulting dispersion relation 
for the effective medium model, as shown in Fig. 
This dispersion relation exhibits classic “avoided cross¬ 
ing” behavior [55] about the resonance frequencies lun 


and ujRH,iso = ■\/7/2 cos- In this model, surface acous¬ 
tic waves (SAWs) in the substrate behave as classical 
Rayleigh waves at frequencies far from the contact res- 

























FIG. 5: Dispersion relation of SAWs in an elastic half space 
coupled to a monolayer of isolated elastic spheres, denoted by 
the blue solid lines. Black dashed lines denote the contact 
resonances, and black dash-dotted lines denote the transverse 
and Rayleigh waves speeds of the substrate. 


onances, and the dispersion curves follow the line corre¬ 
sponding to the substrate Rayleigh wave speed cr m- 
Conversely, sphere motions dominate those of the sub¬ 
strate at frequencies close to the contact resonances. For 
phase velocities greater than cy, which correspond to the 
region w > ct^, the wave numbers that solve Eq. (13) 
are complex valued; these solutions are “leaky” modes 
for which energy leaves the surface of the substrate, and 
radiates into the bulk. This isolated spheres case is par¬ 
ticularly applicable in systems where adhesion between 
particles is negligible, e.g. for: macroscale particles with¬ 
out lateral compression where the dominant static com¬ 
pression is due to gravity and is between the particles 
and substrate; or for microscale particles, if the separa¬ 
tion distance between particles is larger than the range 
of adhesion forces. 


2. Interacting Spheres 



FIG. 6: (a) Blue solid lines denote the dispersion relation of 
SAWs in an elastic halfspace conpled to a monolayer of inter¬ 
acting elastic spheres. Black dashed lines denote the contact 
resonances, and black dash-dotted lines denote wave speeds in 
the substrate, (b)-(e) relative amplitudes of the displacement 
variables uq (black solid lines), Wo (red solid lines), Q (black 
dotted lines), Z (red dotted lines), and R9 (blue dotted lines), 
corresponding to the branch denoted by the same numeral in 
(a). The amplitudes are normalized such that the sum of the 
squares is unity. 


In Fig. [^a), we plot numerical solutions of Eq. ( [I4| ) for 
the long wavelength limit of the effective medium model 
with interacting spheres. The amplitudes of the sphere 
and substrate displacements are calculated in the same 
manner as in Fig. and are plotted in Fig. I^b-e). In 
Fig. I^a), we observe features qualitatively similar to the 
dispersion relation for isolated spheres in Fig. with the 
exception of a third avoided crossing at frequency lohr- 
The mode shapes reveal the ways in which each of the 
branches are influenced by the contact resonances, as well 
as long and short wavelength asymptotic behavior of our 
system. In the long wavelength limit, the substrate mo¬ 
tions closely resemble Rayleigh SAWs with a mix 


of vertical and horizontal motions. Since the frequen¬ 
cies of waves in this regime are well below the contact 
resonances, the effect of the spheres is negligible, and 
the monolayer moves in phase with the substrate sur¬ 
face. At short wavelengths, it can be clearly seen that 
the first, second, and third lowest branches exhibit mo¬ 
tions dominated by the displacements Q, Z, and 9, re¬ 
spectively (each corresponding to a resonant mode of the 
monolayer), while the highest branch tends toward the 
Rayleigh SAW. The effects of proximity to the contact 
resonances are well illustrated, for example, by branch III 
of Fig. [^a), which exhibits large vertical sphere motions 
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at its starting point near wat, resembles the Rayleigh 
SAW as it approaches and crosses the cr line, and tran¬ 
sitions into large rotational sphere motions after bending 
around the avoided crossing with ujrh- 

In order to examine the behavior of our system 
throughout the entire Brillouin zone, we superimpose 
the dispersion curves for the effective medium model of 
interacting spheres on an elastic base including higher 
order spatial derivative terms (the full Eq. (HI) with 
the dispersion curves for the discrete monolayer on a 
rigid substrate (Eq. as shown in Eig. At long 

wavelengths the discreteness of the monolayer is insignif¬ 
icant, and the dispersion is well described by the effective 
medium model. Eurthermore, we note that at long wave¬ 
lengths the dispersion curves calculated using the effec¬ 
tive medium model including higher order terms shown 
in Fig. only slightly deviates from the dispersion cal¬ 
culated using the effective medium model with the higher 
order terms neglected shown in Fig. I^a). The only no¬ 
ticeable effect is a downshift in frequency of the avoided 
crossing between the Rayleigh wave and the ujrh res¬ 
onance; since the latter intersects at the highest wave 
vector of the three contact resonances, calculations with 
Eq. (141 in this case are the least accurate. In Fig. ^ at 
short wavelengths, beyond the avoided crossings with the 
Rayleigh wave branch, the elasticity of the substrate has 
little effect on the dynamics, and the dispersion can be 
described using the discrete model for interacting spheres 
on a rigid substrate. We suggest that by “stitching to¬ 
gether” the effective medium model for spheres on an 
elastic substrate with the discrete model for spheres on a 
rigid substrate, we can simultaneously capture the inter¬ 
action of SAWs with the monolayer at long wavelengths 
and effects caused by the discreteness of the spheres at 
short wavelengths. Past the avoided crossings, the two 
sets of curves in Fig. stitch together smoothly, result¬ 
ing in a full picture of the monolayer dynamics on the 
elastic substrate. 


V. EXPERIMENTAL IMPLICATIONS 

We expect the presented results to be useful for pre¬ 
dicting complex dynamic responses and extracting effec¬ 
tive contact stiffnesses from measurements of acoustic 
dispersion in a manner similar to Boechler et al. |I2j . 
The findings described above invite several questions, in¬ 
cluding whether our model of a square lattice is appli¬ 
cable to results on hexagonally packed monolayers, and 
why horizontal-rotational resonances were not observed 
in the experiment |12j . 

We believe that the assumption of the square lattice 
is not essential. For isolated spheres, Eqs. Q and Q 
with Gn and Gs set to zero can be obtained for any ar¬ 
rangement of the spheres, periodic or random, with the 
only parameter depending on the arrangement being the 
surface area per sphere A. For interacting spheres, the 




FIG. 7: Blue lines denote the SAW dispersion relation with 
spatial derivative terms inclnded and red lines denote the dis¬ 
persion relation of a discrete monolayer adhered to a rigid 
base. Solid and dotted lines denote, respectively, valid and 
invalid ranges for the two models. Black dashed and dash- 
dotted lines are the same as in Fig. [^a). 


results generally do depend on the lattice structure and 
the propagation direction. However, the contact reso¬ 
nances given by Eq. (12) correspond to the fc = 0 limit 
and, consequently, do not depend on the propagation di¬ 
rection. The relative positions of the three contact res¬ 
onances may be different in the long wavelength limit 
between a hexagonal and square packed lattice, but their 
presence should still be expected in both cases. 


We suggest that the reasons why horizontal-rotational 
resonances were not observed in Ref. |12j may be the 
following. Since the measurements were not sensitive to 
horizontal motion, the lorh and ujrr resonances could 
only be detected when they hybridized with SAWs near 
avoided crossings, and since the avoided crossings with 
UJRH and ujhr resonances are more narrow than the one 
with the ujjq resonance, they could have been missed. 
Furthermore, our model assumes that all spheres are ei¬ 
ther connected by identical springs or are isolated. If 
the contact stiffness between spheres were to vary widely 
(some neighboring spheres being in contact and others 
not, for example), then distinct resonances may be ab¬ 
sent. In addition, the upper {ujrh) resonance may have 
been outside the range of the measurements in Ref. m- 
Further experimental studies of monolayer dynamics in 
conjunction with exploration of ways to control sphere- 
to-sphere contacts should help resolve the discrepancy 
between the theory and experiment. 


While the main focus of this work has been on mi- 
crogranular monolayers, our theory is equally valid for 
macroscale systems. In this case, the contact springs 
would be determined by gravity and, possibly, applied 
lateral static compression [^, rather than by adhesion 
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forces. We note that several past experimental works on 
macroscale granular systems [9] have observed system¬ 
atic upshifts in frequency relative to theoretical predic¬ 
tions, and have suggested uncertainties in material pa¬ 
rameters and experimental setups, as well as deviations 
from Hertzian contact behavior as possible causes. As 
per the results from our model, the presence of additional 
degrees of freedom and interactions between spheres and 
substrate may also be possible causes. In the absence of 
the external lateral compression, highly nonlinear “sonic 
vacua” [3] should also be expected. Generally, as ampli¬ 
tudes are increased, interesting nonlinear dynamics are 
expected for both micro- and macroscale monolayers due 
to nonlinearity of Hertzian contacts between the particles 
[am and between the particles and the substrate m- 

VI. CONCLUSION 

We have developed a model for wave propagation in 
granular systems composed of a monolayer of spheres on 
an elastic substrate. Our model expands on those used 
in previous works by including the elasticity of the sub¬ 
strate, horizontal and rotational sphere motions, shear 
coupling between the spheres and substrate, and inter¬ 
actions between adjacent spheres. We have shown that 
a monolayer of interacting spheres on a rigid substrate 
supports three modes involving vertical, horizontal, and 
rotational motion. In the long-wavelength limit, these 
modes yield three contact resonances, one purely ver¬ 
tical and two of mixed horizontal-rotational character. 
On an elastic substrate, these resonances hybridize with 
the Rayleigh surface wave yielding three avoided cross¬ 
ings. For isolated spheres, the frequency of the lower 
horizontal-rotational resonance, in the absence of bend¬ 


ing rigidity, tends to zero and only two contact resonances 
with two respective avoided crossings remain. 

By comparing the effective medium (valid for long 
wavelengths) to the discrete formulation of our model, 
we have demonstrated that for the presented microsphere 
monolayer example, the effective medium model can be 
used to describe the interaction of the contact resonances 
with the Rayleigh waves in the substrate, but loses accu¬ 
racy at shorter wavelengths. In that case, the substrate 
can be considered rigid, and the discrete model is more 
appropriate. This model is scalable in that it can be 
adapted for use with both macro- and microscale sys¬ 
tems, and provides a means to experimentally extract 
contact stiffnesses from dynamic measurements. Oppor¬ 
tunities for future studies include exploration of anal¬ 
ogous models for granular monolayers in the nonlinear 
regime, as well as analysis of the transverse modes of a 
monolayer of spheres on an elastic substrate (the latter 
involves transverse horizontal displacement and rotations 
of the spheres, as well as shear horizontal acoustic waves 
in the substrate). Further experiments with macro- and 
microscale granular monolayers will help guide the mod¬ 
eling effort. 
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